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INTRODUCTION  AND  PROBLEM  STATEMENT 


Some  ultra  high  pressure  vessels  use  threaded  connectors  as  end 
enclosures*  Furthermore,  it  is  sometimes  necessary  to  remove  the  end 
enclosure  quickly,  thus,  the  sectored  thread.  Figure  1,  has  been  developed. 
With  this  design,  the  end  cap  can  be  removed  with  only  a  single  one-quarter 
turn  rather  than  several  turns.  When  it  is  necessary  to  maximize  the  load 
bearing  capability  of  the  connection,  the  thread  area  must  be  maximized  and 
generous  root  radii  at  the  ends  of  the  thread  sectors  are  not  possible.  This 
report  summarizes  the  results  of  a  simulation  fatigue  test  program  to  evaluate 
the  effects  of  small  root  radius  and  surface  treatment  on  the  fatigue  crack 
initiation  life. 

The  vessel  in  consideration  has  four  sectors  and  the  dimension  of  the 
root  radius  is  0.010  in.  R  min  (0.254  mm  R  min)  to  0.020  in.  R  max  (0.508  mm 
R  max).  The  material  was  pressure  vessel  steel  ASTM  A723  Grade  I  Class  4  with 
nominal  yield  strength  of  170  ksi  (1172  MPa).  The  inner  diameter  (ID)  was 
5.165  in.  (1.31  cm),  the  major  outer  diameter  (OD)  (2b^)  was  8.75  in.  (22.2 
cm),  and  the  minor  0D  (2b2)  was  8.41  in.  (21.4  cm).  The  cylinder  was  designed 
to  contain  72  ksi  (496  MPa)  internal  pressure  and  was  swage  autof rettaged  to 
produce  the  100  percent  overstrain  condition. 

Three  different  treatments  of  the  notch  root  radius  were  studied.  First, 
milling  to  the  required  dimensions  with  a  faced  cutting  tool;  second,  milling 
followed  by  electropolishing  to  both  Increase  the  root  radius  somewhat  and  to 
produce  a  smoother  surface  finish;  and  third,  milling  followed  by  glass  ball 
peening.  Steel  shot  could  not  be  used  because  of  the  very  small  root  radius. 
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An  earlier  study  (ref  1)  has  shown  that  electropolishing  a  grooved  sample 
can  substantially  increase  fatigue  life.  The  as-machined  grooves  were  loaded 
in  a  zero-to-tension  manner  such  that  failure  occurred  in  about  20,000  cycles. 
Applying  the  same  loads  to  samples  that  had  been  machined  and  electropolished 
increased  fatigue  life  to  about  35,000  cycles.  This  Increase  was  attributed 
to  the  removal  of  material  damaged  during  machining  by  the  electropolishing 
process.  Although  the  problem  considered  here  is  quite  different  (a  tension- 
tension  loading  and  much  lower  life) ,  electropolishing  is  worthy  of 
Investigation  at  this  time. 


EXPERIMENTAL  PROCEDURE 

Previous  studies  (refs  2,3)  have  shown  that  adequate  simulation  of  OD 
initiated  failures  can  be  obtained  by  using  the  specimen  shown  in  Figure  2. 
The  basic  assumption  made  in  using  this  specimen  is  that  crack  Initiation 
occurs  at  the  notch  due  to  local  stresses  regardless  of  how  the  stresses  are 
produced.  With  the  simulation  specimen,  both  the  tensile  residual  stress  and 
the  alternating  stress  due  to  internal  pressurization  can  be  applied  by 
adjusting  the  load  P. 


1-V.  P.  Greco  and  B.  L.  Pennell,  "Low  Cycle  Notched-Fatigue  Behavior  of 
Broached  and  Electrochemically  Treated  Steels  (Electropolishing  and  Chromium 
Plating),"  presented  at  the  56th  Annual  Convention,  American  Electroplater's 
Society,  Detroit,  MI,  June  1969. 

^J.  A.  Kapp  and  J.  H.  Underwood,  "Service  Simulation  Tests  to  Determine  the 
Fatigue  Life  of  Outside-Diameter-Notched  Thick-Wall  Cylinder,"  Experimental 
Mechanics ,  Vol.  22,  No.  3,  March  1982,  pp.  96-100. 

^J.  H.  Underwood,  "Fatigue  Life  Analysis  and  Tensile  Overload  Effects  With 
High  Strength  Steel  Notched  Specimens,"  High  Pressure  in  Science  and 
TechnoloBV.  Part  TT  :  FlnlHa.  P.ns1np»r{m>.  anH  Safotv.  Vol  .  22 


The  former  investigations  modeled  a  somewhat  different  OD  notch 
configuration  that  was  substantially  deeper  than  the  problem  considered  here. 
In  those  cases,  the  crack  achieved  critical  size  immediately  upon  initiation 
in  both  the  pressure  loading  and  simulation  loading  conditions.  Thus,  the 
fatigue  crack  initiation  life  was  determined  as  the  number  of  cycles-to- 
failure.  The  cylinder  we  studied  is  somewhat  different.  Autofrettage  and 
internal  pressurization  would  result  in  failure  shortly  after  crack 
initiation,  but  the  simulation  specimen  can  accommodate  large  amounts  of 
subcritical  crack  growth.  To  alleviate  this  discrepancy,  a  small  flat  was 
produced  on  the  ID  of  the  simulation  specimen  directly  in  line  with  the  OD 
notch.  An  ultrasonic  probe  was  placed  there  to  monitor  crack  initiation. 
Cycling  proceeded  until  a  crack  about  0.05  in.  (1.3  mm)  deep  was  observed. 
This  depth  was  chosen  because  it  was  the  minimum  depth  that  could  be  reliably 
detected.  Once  the  crack  was  indicated,  the  notch  was  cut  from  the  specimen 
and  broken  apart  to  reveal  the  fracture  surface.  The  depth  of  the  fatigue 
crack  was  then  measured. 

The  fatigue  crack  initiation  life  was  determined  by  calculating  the 
number  of  cycles  necessary  for  the  crack  to  grow  the  measured  amount.  The 
calculated  cycles  were  then  subtracted  from  the  total  measured  cycles.  This 
calculation  was  based  on  the  assumption  that  once  the  crack  initiates,  it  has 
dimensions  equal  to  the  notch  depth  (b^-b2) ,  and  the  simulation  specimen  acts 
as  a  solid  unnotched  ring  with  the  OD  equal  to  2b^.  Stress  intensity  factors 
for  the  smooth  ring  subject  to  the  imposed  loading  were  known  from  Gross  and 
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Srawley  (ref  4),  and  crack  growth  rates  were  determined  in  a  companion  study. 


STRESS  ANALYSIS 

The  basic  prerequisite  of  the  simulation  tests  is  to  produce  the  same 
stress  in  the  root  of  the  notch  by  simulation  loading  as  that  which  occurs  in 
the  autof rettaged  and  pressurized  cylinder.  This  was  accomplished  by  using 
finite  element  stress  analysis  for  one  notch  configuration,  and  a  Neuber 
diagram  (ref  5).  The  finite  element  results  for  the  single  notch 
configuration  were  used  to  determine  the  relative  error  in  the  Neuber  estimate 
for  the  same  configuration,  which  was  applied  to  the  remaining  Neuber  results. 
The  stress  concentration  factors  determined  this  way  were  used  to  estimate 
crack  initiation  life  using  a  "local  strain"  method  discussed  below. 

Three  loading  conditions  were  evaluated  which  required  two  finite  element 
models.  Internal  pressure  and  autof rettage  loading  were  studied  using  the 
finite  element  mesh  in  Figure  3(a)  and  simulation  loading  was  analyzed  with 
the  mesh  shown  in  Figure  3(b).  The  smaller  mesh  was  used  for  pressure  and 
autofrettage  loading  because  of  the  eightfold  symmetry  of  the  cylinder.  There 
are  no  planes  of  symmetry  in  the  simulation  specimen;  the  entire  structure 
must  be  studied.  The  refinement  of  the  meshes  around  the  root  radius  is  shown 
in  Figure  3(c) . 

^B.  Gross  and  J.  E.  Srawley,  "Analysis  of  Radially  Cracked  Ring  Segments 
Subject  to  Forces  and  Couples,"  ASTM  STP  632,  American  Society  for  Testing 
and  Materials,  1977,  pp.  39-56. 

5A.  C.  Ugural  and  S.  K.  Fens  ter,  Advanced  Strength  and  Applied  Elasticity, 
American  Elsevier,  New  York,  1975,  p.  83. 


The  analysis  was  conducted  using  cubic  isoparametric  elements  in  the  APES 
finite  element  program  (ref  6).  Because  these  powerful  elements  were  used,  a 
relatively  coarse  mesh  was  possible.  For  pressure  and  autof rettage ,  166  nodes 
were  used  and  for  the  simulation  specimen,  229  nodes  were  necessary.  It  is 
estimated  that  to  obtain  equivalent  accuracy  using  constant  strain  elements, 
about  4000  degrees  of  freedom  would  be  necessary  for  the  pressure  and 
autofrettage  loading  and  about  5500  degrees  of  freedom  would  be  required  for 
the  simulation  specimen  (ref  6) . 

The  loading  of  these  models  is  straightforward  except  for  the 
autofrettage  case.  The  residual  stresses  were  simulated  using  the  thermal 
loading  analogy  of  Hussain  et  al  (ref  7).  It  was  noted  in  Reference  7  that 
autofrettage  residual  stresses  have  exactly  the  same  distribution  as  if  the 
cylinders  were  subjected  to  a  logarithmically  varying  axisymmetric  temperature 
distribution  in  the  same  smooth  cylinder.  By  loading  a  notched  finite  element 
model  with  the  same  temperature  distribution  as  that  which  produces  100 
percent  overstrain  in  a  smooth  cylinder  with  an  0D  of  2b^,  the  effects  of  the 
notch  on  the  residual  stresses  can  be  accurately  determined  (refs  7,8). 


^1.  N.  Gifford,  Jr.,  "APES  -  Second  Generation  Two-Dimensional  Fracture 
Mechanics  and  Stress  Analysis  by  Finite  Elements,”  Report  4799,  Naval  Ship 
Research  and  Development  Center,  1979. 

^M.  A.  Hussain,  S.  L.  Pu,  J.  D.  Vasilakis,  and  G.  P.  O'Hara,  Journal  of 
Pressure  Vessel  Technology,  Vol.  102,  No.  3,  1980,  p.  314. 

8S .  L. Pu  and  M.  A.  Hussain,  Journal  of  Pressure  Vessel  Technology,  Vol.  103, 
May  1981,  p.  301. 
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LOCAL  STRAIN  LIFE  ESTIMATE 


Accurate  predictions  of  the  fatigue  life  of  notched  components  have  been 
made  using  the  "local  strain"  method  (ref  9).  This  technique  works  as 
follows:  the  strain  range  (Ae)  and  strain  ratio  (R  *  Snin/Smax)  at  the  root 

of  the  notch  are  estimated  using  an  elastic  analysis  of  the  specimen.  Crack 
initiation  life  is  estimated  from  the  steel  strain-life  property  generated  by 
using  strain-controlled  loading,  and  by  accounting  for  the  effects  of  R.  The 
equivalent  local  strain  range  was  calculated  by 

( amin“amax) 

Ae  =  kf -  (1 

E 


and  the  strain  ratio  was 

^  “  °min/ (Jmax  (2) 

Where  kf  is  the  fatigue  stress  concentration  factor  which  is  related  to  the 
static  elastic  stress  concentration  factor  kt  for  this  material  by  (ref  9): 

kf  -  1  +  qOct-l)  (3) 


q  -  — — - 
1+a/  p 


(4) 


_  300 

a  =»  0.001( - ) 


1.8 


(5) 


2078 

a  =>  0.0254( - ) 


1.8 


(6) 


^H.  S.  Reemsnyder,  "Constant  Amplitude  Fatigue  Life  Assessment  Models,"  SAE 
Technical  Paper  Series  820688,  Proceedings  of  the  SAE  Fatigue  Conference, 
April  1982,  pp.  119-132. 


where  o0  is  the  ultimate  tensile  strength  of  the  material  in  units  of  ksi  in 
Eq.  (5)  or  in  units  of  MPa  in  Eq.  (6). 

RESULTS  AND  DISCUSSION  OF  RESULTS 

The  stress  analysis  results  are  summarized  in  Table  I.  The  stress 
concentration  factors  (kt)  determined  from  the  finite  element  models  for  each 
loading  condition  with  the  minimum  root  radius  are  reported.  Also,  kt  for 
several  notch  geometries  within  the  tolerance  are  given  from  a  Neuber  diagram. 
As  can  be  seen  the  agreement  is  quite  good,  but  the  finite  element  results  are 
somewhat  lower  than  that  predicted  by  the  Neuber  diagram.  Since  the  Neuber 
diagram  gives  an  estimated  kt  for  a  geometry  that  is  somewhat  different  than 
the  actual  cylinder,  it  is  assumed  that  the  finite  element  results  are  more 
accurate  and  should  be  used  in  subsequent  analyses.  As  kt  is  not  a  constant 
with  respect  to  loading  condition,  an  average  was  determined  from  the  finite 
element  solutions.  This  appears  in  the  final  column  of  the  Table  I. 

To  estimate  kt  for  the  more  generous  root  radii,  it  was  decided  to  use 
only  the  Neuber  diagram.  The  finite  element  solutions  were  too  time 
consuming.  To  improve  the  accuracy  of  the  Neuber  diagram,  a  ratio  was 
established  from  the  one  condition  in  which  both  finite  element  and  Neuber 
results  were  available.  This  ratio  was  then  multiplied  by  the  kt  estimated 
from  the  Neuber  diagram.  These  values  of  kt  appear  as  the  final  column  in  the 
table  as  a  function  of  root  radius. 

The  kt  values  determined  in  this  manner  were  used  to  estimate  the  life  of 
the  notched  cylinder  using  the  local  strain  concepts.  Equations  (4)  and  (5) 
were  used  along  with  other  information  given  in  Table  II  to  generate  first  kf 
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Chen  Che  esCimaCed  local  sCraln  range.  These  are  also  given  in  Table  II. 
Knowing  boch  Che  sCrain  range  and  Che  scrain  racio,  an  escimace  of  Che  number 
of  cycles  Co  iniCiaCe  a  crack  can  be  made  wich  Che  help  of  Figure  4.  The  e-N 
curve  for  Che  same  maCerial  CesCed  here  was  generaCed  previously  in  a 
companion  scudy,  as  was  Che  crack  growch  race  as  a  funcclon  of  scress 
inCensiCy  facCor  range  and  R  (Figure  5).  The  number  of  cycles  necessary  Co 
IniCiaCe  a  crack  are  given  in  Che  final  column  of  Table  II.  These  resulcs 
seem  as  expecCed,  namely  as  Che  rooC  radius  increases,  Che  local  loading  is 
less  severe  and  longer  lives  should  be  observed. 

The  measured  facigue  lives  of  Che  specimens  CesCed  are  summarized  in 
Table  III.  In  chis  Cable  all  of  Che  informaCion  gaChered  on  Che  simulaCion 
specimen  Is  given:  Che  measured  rooc-radius  ( p) ,  Che  measured  lengch  of  Che 
facigue  crack  chac  developed  (Ac),  Che  CoCal  number  of  cycles  accumulaCed  when 
Che  Cesc  was  Cerminaced  (Nt) ,  Che  calculaCed  number  of  cycles  Co  propagace  a 
crack  co  chaC  depCh  (N^c) ,  and  Che  esCimaCed  number  of  cycles  Co  iniCiaCion 
(N^  »  Nt  -  Nac).  Nac  was  deCermlned  wich  Che  help  of  Figure  5  and  Che  sCress 
inCensiCy  facCors  for  a  curved  bar  wich  an  exCernal  crack  subjecC  Co  combined 
bending  and  compression  (ref  4).  As  sCaCed  above,  ic  was  assumed  chac  ac 
IniCiaCion,  che  simulaCion  specimen  acCs  as  a  smooch  curved  specimen  wich  an 
exCernal  crack.  Ics  dimensions  were  assumed  Co  be:  OD  =■  2b^  *  8.75  in. 

(22.2  cm),  ID  =  5.17  in.  (13.12  cm),  and  crack  lengch  c  *  b^  -b2  *  0.17  In. 


+B .  Gross  and  J.  E.  Srawley,  "Analysis  of  Radially  Cracked  Ring  SegmenCs 
SubjecC  Co  Forces  and  Couples,"  ASTM  STP  632,  American  SocieCy  for  TesClng 
and  MaCerials,  1977,  pp.  39-56. 


(0.43  cm).  From  Reference  4,  it  was  found  that  the  stress  intensity  factor 
range  (AK)  for  this  initial  condition  was  about  49  ksi/in.  (54  MPa/m).  After 
as  much  as  0.070  in.  (0.18  cm)  of  crack  extension,  AK  was  about  52  ksi/in.  (57 
MPa/m).  Since  AK  changed  little,  small  errors  were  incurred  when  the  cycles 
to  extend  the  crack  were  calculated  as  N^.  -  Ac/(dc/dN(AK  *  50.5  ksi/in.)) , 
(55.5  MPa/m) •  From  Figure  5  at  AK  ■  50.5  ksi/in.  (55.5  MPa^m)  and  R  *  +0.5, 
dc/dN  *  2.1  x  10-i*  in. /cycle  (5.3  x  10“6  m/cycle).  The  calculated  values  of 
NAc  were  small  and  probably  very  approximate  because  of  the  assumptions  made. 
If  more  were  known  of  the  growth  of  small  fatigue  cracks  near  a  notch,  a 
better  estimate  of  NAc  could  be  made.  Although  the  values  reported  in  the 
table  are  only  estimates,  they  are  still  reported  to  account  for  the  fact  that 
the  fracture  surfaces  showed  fatigue  cracks  of  different  lengths  and  certainly 
some  number  of  cycles  were  necessary  to  grow  the  cracks. 

The  data  are  further  compared  graphically  in  Figure  6.  From  this  figure 
we  see  that  the  local  strain  life  estimate  gives  a  lower  bound  of  the  measured 
behavior,  especially  with  the  larger  root  radii.  The  model  can  underestimate 
life  by  as  much  as  50  percent  in  the  as-machined  samples.  This  difference 
can  be  attributed  to  scatter.  The  important  observation  to  make  from  this 
plot  is  that  changing  the  root  radius  will  affect  fatigue  crack  initiation 
life  more  than  surface  treatment.  For  example,  local  strain  analysis  suggests 
that  the  life  to  initiate  a  crack  at  a  notch  with  a  root  radius  of  0.020  in. 

^B.  Gross  and  J.  E.  Srawley,  "Analysis  of  Radially  Cracked  Ring  Segments 
Subject  to  Forces  and  Couples,"  ASTM  STP  632,  American  Society  for  Testing 
and  Materials,  1977,  pp.  39-56. 


(0.508  urn)  should  be  about  twice  as  along  as  to  initiate  a  crack  with  a  root 
radius  of  0.010  in.  (0.254  mm).  The  experimental  data  shows  that  a  somewhat 
greater  ratio  is  obtained.  Comparing  specimens  18  and  9  with  root  radii  of 
0.012  in.  (0.305  mm)  and  0.018  in.  (0.457  mm)  respectively,  life  was  increased 
by  more  than  a  factor  of  three.  The  data  from  the  electropolished  samples 
cannot  be  distinguished  from  the  as-machined  samples,  nor  can  the  glass 
peenlng  results  be  interpreted  as  improving  fatigue  life.  This  relative 
insensitivity  of  fatigue  life  to  notch  surface  condition  is  in  agreement  with 
prior  results  (ref  3),  of  fatigue  life  tests  with  notch  root  stresses  greatly 
in  excess  of  the  yield  strength.  There  was  one  electropolished  sample 
(specimen  no.  2)  that  indeed  indicated  an  improvement,  but  since  that  result 
could  not  be  repeated,  it  was  assumed  to  be  an  artifact  and  inconclusive. 

Although  the  results  of  changing  surface  treatments  were  ineffective  in 
this  case,  we  cannot  generalize  these  results.  We  tested  under  only  one 
loading  condition  resulting  in  a  very  low  life.  It  is  clear  that  under  the 
conditions  studied,  any  effect  due  to  either  electropolishing  or  glass  ball 
peenlng  were  overcome  by  the  substantial  inelastic  behavior  at  the  root 
radius.  It  is  well  known  that  peening  increases  fatigue  life  in  the  high 
cycle  regimen  and  previous  studies  (ref  1)  also  indicate  the  same  for 
electropolishing . 

*-V.  P.  Greco  and  B.  L.  Pennell,  "Low  Cycle  Notched-Fatigue  Behavior  of 
Broached  and  Electrochemically  Treated  Steels  (Electropolishing  and  Chromium 
Plating),”  presented  at  the  56th  Annual  Convention,  American  Electroplater's 
Society,  Detroit,  MI,  June  1969. 

^J.  H.  Underwood,  "Fatigue  Life  Analysis  and  Tensile  Overload  Effects  With 
High  Strength  Steel  Notched  Specimens,"  High  Pressure  in  Science  and 
Technology,  Part  II;  Fluids,  Engineering,  and  Safety,  Vol.  22,  1984,  pp. 
209-214. 
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CONCLUSIONS 


The  low  cycle  fatigue  crack  initiation  life  of  an  autof rettaged  cylinder 
with  a  notch  on  its  OD  can  be  adequately  determined  with  a  simulation  specimen 
if  care  is  taken  to  detect  crack  initiation.  The  effects  of  the  root  radius 
of  the  notch  can  be  easily  determined  experimentally  with  this  technique.  The 
result  is  that  increasing  the  root  radius  from  about  0.010  in.  (0.254  mm)  to 
about  0.020  in.  (0.508  mm)  increases  the  initiation  life  by  about  a  factor  of 
three.  Whether  the  cylinder  is  cycled  as-machined,  or  followed  by  electro¬ 
polishing  or  glass  ball  peening,  there  is  no  effect  on  fatigue  crack 
initiation  life.  Further  testing  with  higher  life  cylinders  should  be 
performed  to  assess  the  effects  of  surface  treatments  more  thoroughly. 

Finally,  the  lower  bound  of  fatigue  initiation  life  can  be  determined  using 
local  strain  analysis,  if  smooth  specimen  strain  range  versus  cycles-to- 
failure  data  are  available. 
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STRESS  ANALYSIS  RESULTS 


54  cm 
89  MPa 


TABLE  II.  ESTIMATED  LIFE  USING  LOCAL  STRAIN  ANALYSIS 


190  ksi  (1310  MPa)  ,  <W  -  ouin 

E  -  30,800  ksi  (212.3  GPa) 


87.1  ksi  (600  MPa), 


r 

i 

i 

i 

Root  Radius 

(in.) 

“i - r 

1  1 

I  ^  1 

q(Eq.  4) 

T - 

1 

j  fcf 

Ae(Eq.  2) 

"1 

1  %  (Cyc! 

1  (Fig.  4 

1 

T 

1 

0.010 

1  1 

1  5.37  I 

I  I 

0.81 

1 

1  4.54 
| 

0.0128 

1 

|  1292 

I 

1 

1 

i 

0.012 

1  1 

1  4.95  1 

1  1 

0.84 

1  4.32 
[ 

0.0122 

I  1756 

1 

i 

1 

i 

0.016 

1  l 

1  4.38  I 

l  1 

0.87 

1  3.94 

0.0111 

1 

I  2268 

i 

1 

| 

0.018 

!  1 

1  4.11  | 

1  1 

0.88 

1  3.73 

0.0105 

I  2387 

1 

i 

1 

0.020 

1  1 

1  3.96  I 

0.89 

1  3.63 

0.0103 

I  2644 

1 

1  -  1 _ 

1 

1 

TABLE  III.  MEASURED  LIVES  USING  THE  SIMULATION  SPECIMEN 


AK  -  50.5  ksi/in.  (55.5  MPa/m)  (Average),  R  »  0.46 
dc/dN  ■  2.1x10"*  in. /cycle  (5.3  x  10“6  m/cycle)  (Figure  5) 
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Figure  6.  Measured  crack  initiation  life  as  a  function  of  notch  root  radius 
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